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Microlensing studies are usually based on the lens equation that is valid only to the 
first order in the gravitational constant G and lens mass M. We consider corrections to the 
conventional lens equation in terms of differentiable functions, so that they can express not 
only the second-order effects of GM in general relativity but also modified gravity theories. 
As a generalization of Ebina et al. (Prog. Theor. Phys. 104 (2000) 1317), we show that, 
provided that the spacetime is static, spherically symmetric and asymptotically flat, the total 
amplification by microlensing remains unchanged at the linear order of the correction to the 
deflection angle, if and only if the correction takes a particular form as the inverse square 
of the impact parameter, whereas the magnification factor for each image is corrected. It 
is concluded that the light curve shape by microlensing is inevitably changed and will thus 
allow us to probe modified gravity, unless a modification to the deflection angle takes the 
particular form. No systematic deviation in microlensing observations has been reported. 
For instance, therefore, the Yukawa-type correction is constrained as the characteristic length 
> 10 14 m. 

§1. Introduction 

Since the pioneering work by PaczynskP and the subsequent detections of MA- 
CHOs,®'® microlensing has been one of the most vital areas in astrophysics. Now, 
it plays an important role also in searching extra-solar planets.®'® Light rays are 
influenced by a local curvature of a spacetime.® It seems natural that the bending 
angle of light rays in modified gravity theories is different from that in general rela- 
tivity. Does such a modified bending angle really makes a change in microlensing? 
We shall examine this problem in this paper. 

Conventional studies of microlensing are usually based on the lens equation that 
maps the source direction into the position on the lens plane, namely the direction 
of each image due to the lens effect.® In particular, the deflection angle is written 
as AGM/bc 2 , where G is the gravitational constant, M is the lens mass, b is the 
impact parameter of the light ray, and c is the speed of light. For the derivation 
of the deflection angle and consequently the lens equation, the post-Minkowskian 
approximation 0(G) or the linear-order approximation O(M) of the Schwarzschild 
metric is employed. 

Do second-order relativistic corrections affect microlensing through the amplifi- 
cation factor? This has been already answered.® The total amplification remains 
unchanged at 0(G 2 M 2 ), whereas the magnification factor for each image is corrected 
at this order. 

The main purpose of this paper is to generalize their relativistic result. In 
particular, theoretical models beyond the theory of general relativity have attracted 
a lot of interests last decades, mostly motivated by the dark energy and dark matter 
problems. In this letter, therefore, we shall reexamine the amplification of lensed 
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images by taking account of corrections in a rather general form, where we assume 
static, spherically symmetric spacetimes. 

We use the units of c = 1 but keep G in order to make iterative calculations 
clearer. 



§2. Corrections to the lens equation and magnification 



2.1. Standard form of lens equation and magnification 

Let us begin by briefly summarizing the derivation of the total amplification for 
microlensingJ^ We denote distances from the observer to the lens, from the observer 
to the source, and from the lens to the source as D^, D$ and D^s, respectively. The 
source and image position angles are denoted by 9 s and respectively. The lens 
equation is written as 

S = 0i- ^a, (1) 

where we used the thin-lens approximation. This gives us a mapping between the 
lens and source planes. We can choose 6$ > without loss of generality. For a point 
mass lens (generally a spherically symmetric lens), the deflection angle a at 0(G) 
becomes 

AGM 

<*=— (2 » 

Hence, the lens equation is 

0s = 0i- |s (3) 
where Be is the angular radius of the Einstein ring as 



9e -'I^Us-- (4) 



For microlensing in our galaxy, the radius is typically 



□ / M \i/2/10kpc\i/2 



where we assumed Dl ~ -Dls- 

The lens equation is solved easily as 

1 



0± = -(05±V05 + 4 ^)- (6) 



We obtain the amplification due to gravitational lensing as 

9! d9j 



A 



9s d9 s 



(7) 



which is a function of the image position 
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For each image at 0±, the amplification factor becomes 

Ai = -^=±\, (8) 
2uVu 2 + 4 2 

where u denotes 9s /9e, the source position in units of the Einstein ring radius. 
For microlensing events, the angular separation between the images is too small to 
measure. All we can measure is the total amplification that is expressed as 

u 2 + 2 

Aotal = - 4 + + - 4 -= r^-j - ( 9 ) 
uy u + 4 

2.2. Modified bending angle and its effects on magnification 

Let us take account of some modification in the bending angle of light rays. 
We focus on modified gravity theories that consider spacetimes to be differentiable 
manifolds. This implies that the modified bending angle can be expressed in terms of 
differentiable functions. We assume also that the lens object produces a static and 
spherically symmetric spacetime, for which the bending angle may take a general 
form as 

AGM . 

a = —[l + F(b)}, (10) 

with a differentiable function F(b) denoting modifications. This point is contrast to 
previous worksP"® in which particular forms such as the second-order of mass in 
the Schwarzschild spacetime and the PPN formalism are assumed. 

For specific models, F(b) may be approximated in terms of power functions bPW^ 
For instance, p = 1/2 for the DGP model as one of brane scenarios,^ and p = 3/2 
for a candidate of massive graviton theories P^'^ 

The lens equation to be solved becomes 

§ s = I -a(0 I ), (11) 
with the modified bending angle as 

a(e I ) = l(i + Ef(\e I \)) : (12) 

where we introduce a nondimensional small quantity e as the expansion parameter 
and f(\9i\) is an arbitrary differentiable function corresponding to F(b). Here, 6$ 
and 9 1 are normalized by 6e as 9 s = 0s/0e = u and 9j = 9j/9e, respectively. We 
should note that /(#/ = +oo) can be absorbed into the mass that is defined at the 
spatial infinity and hence f(9j = +oo) can be taken as zero. Henceforth, we assume 
/(6»i) — y as Oj — y oo. 

This lens equation is rewritten as 

9j-9 s 0i-l-ef(\9j\) = O. (13) 

The lens equation (|13|) can be solved iteratively in terms of e. For 9j > 0, the 
lens equation is rewritten as 

9j - 6 s h ~ 1 - = 0. (14) 
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A perturbative form as 

9' + = e + + ec/) + + 0(e 2 ), (15) 

is substituted into the image position 6j in Eq. (|14p . where the prime denotes a 
quantity including effects by the modified bending angle. Here, one can assume 
9 + > because of > 0, within the limit that perturbative calculations hold. 
Then, we find the solution as 

' ~ +e _^+L + ( £ 2 ). (16) 
9f + 4 

For 6j < 0, the lens equation becomes 

Bj - h0 1 - 1 - ef(-h) = 0. (17) 
In the similar manner to the positive case, one finds the solution as 



'1 + 4 



+ 0(0- (18) 



By using Eqs. (|16p and (|18p . we obtain the amplification for each image as 

u 2 + 2 , (\ „ \ ^ /J2 



A ^^Tt ± [^) +0(s) - (19) 

where we define the correction term £± as 

4,4*^f^ + ^) T /, ± j t )^. ,20) 

Here, we assume that the sign of the Jacobian for the lens mapping does not change 
at the linear order of the modified bending angle as far as perturbative treatments 
are valid. Clearly the modified bending angle leads to a change in the amplification 
for each image. 

By using Eq. (|19|) . we obtain the total amplification as 

A[ otal = + e(£ + - tJ) + 0(e 2 ). (21) 

According to Eq. ([2ip . l + = £_ is equivalent to the condition that the total 
amplification for the modified bending angle remains unchanged at the linear order. 
By using Eq. (|2U|) . this condition is rewritten as 

d e + f(e+) _ d (-6L)/(-0_) 



du Vn 2 + 4 du ^J u 2 + 4 
where we should note that 9± is a function of u. This is integrated as 

s + f(S + ) + U(-o-) =ccmst (23) 
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For the asymptotically flat case, this constant must vanish as shown below. In 
the limit as u — > oo, we have 9+ — > u — > oo and —> 0. Hence the L.H.S. of Eq. 
(|23|) becomes /(#+ = oo) = as u — > oo, which means that the above constant 
vanishes. Therefore, the condition of 1+ = i- is expressed as 

§+/(§+) = -Lf (-$-). (24) 

By using the identity as 6+9- = —1, this is rewritten simply as 

9^)=9{^j, (25) 

where we define g{x) = xf(x). Let us prove that g(x) is a constant. We expand the 
differentiable function g{x) in the Laurent series as 

oo 
r=— oo 

where a r is some constant and singular points may exist. This gives also g(l/x) as 
For these expansions, Eq. ([25]) tells us 



r=— oo 



a r = (i-r, (28) 
which allows us to rearrange the expansion of g{x) as 

g(x) =a + a 1 (^x + ^j +a 2 (^x 2 + +■■■ . (29) 

The asymptotic flatness requires that the bending angle vanishes at the spatial 
infinity. Namely, x~ 1 f(x) vanishes as x — > oo, which means a r = for r > 2. Then, 
g{x) = ao + a\{x + x" 1 ) gives the bending angle defined by Eq. (fl~2j) as 



5(6) = i 
b 



l + ^T + 01 ( 1 + P 



(30) 



However, the two terms with a± are rewritten as 



6»i 4GMD L D LS , 

1 + z>,^ " (31) 

This expression makes no sense, because the bending angle is calculated by assuming 
the observer and source located at the null infinity and hence it includes neither D^s, 
Dl nor Ds- Therefore, a\ must vanish. It should be noted that this subtle argument 
is not true of Eq. (|26p , since there are a lot of ways for combining terms in an infinite 
series. 
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As a result, we find only g(x) = ao, which leads to f(x) = clq/x. This equation 
states exactly that the modified bending angle is proportional to the inverse square 
of the impact parameter (e.g., in Eq. (|12[) ). One example of such inverse-square 
corrections is the second-order approximation of the Schwarzschild lens.^'^ An- 
other example in general relativity is the Reissner-Nordstrom solution representing 
a charged black hole, for which the correction to the bending angle is ~ q 2 jb 2 for a 
small charge q. 

Note that there is a subtle but large difference between treatments of a\ and ao- 
This is because the deflection angle depends on the lens mass and impact parameter 
but not any other distances such as Dls in the lensing theory based on the null 
infinity condition.^ The coefficient of a\ is given by Eq. as the ill-defined sum 
of the (nondimensional) unity and the combination of the mass and distances, so 
that the distances cannot be removed by using the impact parameter. On the other 
hand, ao has a coefficient as the inverse of the impact parameter squared. Hence, 
the coefficient is well-defined and thus allowed. 

2.3. Discussion 

Do modified magnification factors which are discussed above have any implica- 
tions for observations? In microlensing, u has a definite dependence on time t, which 
is expressed as u = {uQ + v^t 2 ) 1 / 2 for uq corresponding to the minimum of the impact 
parameter and vt denoting the transverse angular velocity of the source with respect 
to the lens. The modified total amplification thus has a different dependence on time 
and thus produces a different shape of the light curve. It should be noted that the 
peak height of the light curve is less informative on modifications, because it depends 
on many parameters and hence there is a degeneracy in parameter spaces. Shapes 
of light curves are of greater importance. In principle, accurate measurements of 
light-curve shapes can be used for detecting (or observationally constraining) the 
modifications discussed above. 

For instance, we take a correction to the deflection angle as the Yukawa-type 
function as 

where K is roughly the inverse of the characteristic length in a modified gravity 
theory or the inverse of a new mass scale. Figures Q] and [2] show effects of the 
modified bending angle on light curves. Figure [1] gives shapes of light curves for two 
cases. One is based on the amplification factor in the standard form and the other is 
for the Yukawa-type modification as f(b) = exp(-Kb). In order to make it clearer 
the deviation from the standard shape of the light curve, we define the ratio as 

SA(u) = ^total I ^total ^ ^ 
Aotal 

where ^total an d ^ ^ a } denote the standard form of total amplification and the 
modified one, respectively. Figure [2] plots 8 A, the deviation from the standard shape 
of the light curve. For these figures, we choose uq = 1, e = 0.2 and K = 1 so that 
effects by the modification can be distinguished by eye. 
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Fig. 1. Light curves for uo = 1: The solid curve is based on the total amplification in the standard 
form. The dashed one includes the Yukawa-type correction with e = 0.2 and K = 1. Time t is 
in units of the Einstein cross time as 9e/vt- 



For this model, a is effectively increased and hence ^total seems to be enhanced. 
However, ^4^ ^ a j is decreased as shown by Fig. [2j This is because the amplification 
is caused by not the bending angle but its derivative. In fact, the Yukawa-type 
correction gives df(b)/db = —Kexp(-Kb) < 0, namely a minus contribution to the 
amplification factor. So far, no systematic deviation in microlensing observations 
has been reported. One can thus put a constraint on K as K~ x > 1000 AU ~ 10 14 
m for e ~ 1, where we assume the photometric accuracy comparable to 0.1 percents 
and b ~ Dl®e ~ O(IAU). This bound is consistent with gravitational experiments 
in the solar system. 

Finally, we mention an asymptotically massless spacetime, for which the above 
theorem still stands by taking M — > and considering the linear order in e. One ex- 
ample is Ellis' wormhole that makes a difference in light curvesp^' since the bending 
angle for this spacetime is ~ b~ 2 at the leading order but ~ b~ 4 at the next order. 



§3. Conclusion 



We investigated corrections to the conventional lens equation in terms of differ- 
entiable functions, so that they can express not only the second-order effects of GM 
in general relativity but also modified gravity theories. It was shown that, provided 
that the spacetime is static, spherically symmetric and asymptotically flat, the total 
amplification by microlensing remains unchanged at the linear order of the correc- 
tion to the deflection angle, if and only if the correction takes a particular form as 
the inverse square of the impact parameter, whereas the magnification factor for 
each image is corrected. It is concluded that the light curve shape by microlens- 
ing is inevitably changed and will thus allow us to probe modified gravity, unless a 
modification to the deflection angle takes the particular form. 

It is left as a future work to use the present formulation to probe modified 
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Fig. 2. The difference of the two light curves in Figure [T] This shows the dependence of 5 A on 
time (through the impact parameter). 



gravity models by microlensing observations. Microlensing in our galaxy is sensitive 
to gravity at short scale around a few AU, whereas cosmological microlensing is more 
useful for the large distance physicsP 
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